Modeling electrolytically top gated graphene 
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We investigate doping of a single-layer graphene in the presence of electrolytic top gating. The 
interfacial phenomena is modeled using a modified Poisson-Boltzmann equation for an aqueous 
solution of simple salt. We demonstrate both the sensitivity of graphene's doping levels to the salt 
concentration and the importance of quantum capacitance that arises due to the smallness of the 
Debye screening length in the electrolyte. 
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I. INTRODUCTION 

Carbon nano-structures show great promise in many 
applications, including chemical and biological sensors. 
While carbon nanotubes (CNTs) have been extensively 
studied in that context for quite some time [3, Q , investi- 
gations of graphene as a sensor are only beginning to ap- 
pear [1, U, Q . Sensory function of carbon nano-structures 
is generally implemented in the configuration of a field 
effect transistor (FET), with a prominent role played 
by the gate potential that controls the current through 
the device. Biochemical applications require good un- 
derstanding of the interaction of carbon nano-structures 
with aqueous solutions [5j, often in the context of the 
electrochemical top gating [6[ . While significant progress 
has been achieved in understanding the interaction of 
CNT-FETs with the electrolytic environment @, H, H, 
similar studies involving graphene have appeared only 
very recently focusing on the screening effect of an 
ion solution on charge transport through graphene based 
FETs *10] , as well as on the measurement of the quantum 
capacitance of graphene as an ultimately thin electrode 
in an aqueous solution [Tlj |. 

The top gating of a graphene based FET with a solid or 
liquid electrolyte presents several advantages compared 
to the conventional back gating with a metallic electrode. 
Upon application of gate voltage, free ions in the elec- 
trolyte re-distribute themselves, forming an electrostatic 
double layer (EDL) at the interface between graphene 
and the electrolytic solution [l2|. Depending on the ion 
concentration, the EDL can be only a few nanometers 
thick, while still providing efficient shielding of graphene. 
As a consequence, the capacitance of the EDL in an elec- 
trolyte can be much higher than the capacitance of the 
back gate, which is typically separated from graphene 
with a layer of Si02 a few hundred nanometers thick [121 ] . 
This property of the EDL enables a much better control 
of the surface potential on the graphene layer, while re- 
quiring a much lower operating voltage that needs to be 
applied to the reference electrode in the electrolyte than 
voltages currently used with back gates [l2j ■ The applied 
voltage then modifies the chemical potential of graphene, 
resulting in changes in its observable properties such as 
conductance. Since properties of the EDL depend on the 



ion concentration, monitoring the resulting changes in 
graphene's conductance can provide a means for sensor 
application, e.g., in measuring the amount of salt in the 
solution. 

On the other hand, referring to the electrical model of 
the electrolytic gating as a series connection of capacitors 
[T3 |. the high gate capacitance in the electrolyte gives a 
much more prominent role to the quantum capacitance 
of graphene than does the back gate [HI [H, [ID, [l6[ . In 
addition, doping levels of an electrolytically top gated 
graphene have been reported recently [6[ to be much 
higher than those obtained with the conventional back 
gate [l7|. At the same time, mobile ions in the solu- 
tion seem to provide a much more effective screening 
of charged impurities underneath the graphene, thereby 
significantly increasing the charge carrier mobility in 
graphene in comparison to some other high-K dielectric 
environments All these facts indicate that elec- 

trolytic top gating provides a means to develop high per- 
formance FETs. 

While the above few experimental observations reveal 
quite fascinating aspects of the graphene-electrolyte in- 
teraction, theoretical modeling of this system seems to 
be lagging behind the experiment. It is therefore desir- 
able, and tempting to discuss doping of a single layer of 
graphene by a remote gate electrode immersed in a thick 
layer of electrolyte by using two simple models: one de- 
scribing graphene's it electron band structure in the lin- 
ear energy dispersion approximation [12j ]. and the other 
describing the distribution of ions in the electrolyte by 
a one-dimensional (ID) Poisson-Boltzmann (PB) model, 
which takes advantage of the planar symmetry of the 
problem [l3j]. However, it should be emphasized that the 
experiments involving electrolytic top-gating of both car- 
bon nanotubes and graphene [ally] use rather hight 
voltages, on the order of 1-2 V, which can cause signifi- 
cant crowding of counter-ions at the electrolyte-graphene 
interface. It is therefore necessary to go beyond the stan- 
dard PB model by taking into account the steric effects, 
i.e., the effects of finite size of ions in the solution. To 
that aim, we shall use the modified PB (mPB) model 
developed by Borukhov et al. [19|, l2fj , which retains ana- 
lytical tractability of the original 1D-PB model. In addi- 
tion, applied voltages beyond 1 V also require taking into 
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account non-linearity of graphene's band energy disper- 
sion, giving small but noticeable corrections to the linear 
approximation. 

We shall consider here a simple 1:1 electrolyte rep- 
resenting an aqueous solution of NaF because both the 
Na + and F~ ions are chemically inert allowing us to ne- 
glect their specific adsorption on the graphene surface 
fill Il3j . In particular, we shall analyze the density of 
doped charge carriers in graphene at room temperature 
(RT) as a function of both the applied voltage and the 
salt concentration to elucidate graphene's sensor ability. 
In addition, we shall evaluate the contributions of both 
graphene and the EDL in the total gate capacitance in 
terms of the applied voltage to reveal the significance of 
quantum capacitance, as well as to elucidate the behav- 
ior of the EDL under high voltages. We shall cover broad 
ranges of both the salt concentration, going from pM to 
a physiologically relevant value, and the applied voltage, 
going up to about 2 V. 

After outlining in section II our theoretical models for 
graphene and the EDL layer, we shall introduce several 
reduced quantities of relevance for these two vastly differ- 
ent systems and present our results in section III. Con- 
cluding remarks are given in section IV. Note that we 
shall use gaussian units (Aireo = 1) throughout the pa- 
per, unless otherwise explicitly stated. 



II. THEORETICAL MODEL 

Graphene is a semi-metal, or a zero-gap semiconduc- 
tor because its conducting and valence ir electron bands 
touch each other only at two isolated points in its two- 
dimensional (2D) Brillouin zone [l§]. The conical shape 
of these bands in the vicinity of these points gives rise 
to an approximately linear density of states, Pl{z) — 
gd\£\/[2TT(?iVF) 2 ], where gd = 4 is the spin and the band 
valley degeneracy factor, and vf ~ c/300 is the Fermi 
speed of graphene, with c being the speed of light in 
vacuum |l2j. In the intrinsic, or undoped graphene, the 
Fermi energy level sits precisely at the neutrality point, 
Ep = 0, also called the Dirac point. Therefore, the elec- 
trical conductivity of graphene is easily controlled, e.g., 
by applying a gate voltage Va that will cause doping of 
graphene's 7r bands with electrons or holes (depending 
on the sign of Va), which can attain the number density 
per unit area, n, with a typical range of n ~ 10 11 to 10 13 
cm~ 2 [12]. In a doped graphene Fermi level moves to 
£f = TiVF^Jit\n\ sgn(n), where sgn(n) = ±1 for electron 
(hole) doping. At a finite temperature T, one can express 
the charge carrier density in a doped graphene in terms 
of its chemical potential p as [2l| 



n(p) = / de p{e) 
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where (3 = (kgT) 1 with kg being the Boltzmann con- 
stant. We shall use in our calculations a full, non- 



linear expression for the n electron band density, p{e), 
given in Eq. (14) of Ref . [l2| . However, for the sake of 
transparency, the theoretical model for graphene will be 
outlined below within the linear density approximation, 
p(e) ~ Pl(s)- We note that this approximation is accu- 
rate enough for low to moderate doping levels, such that, 
e.g., \p\ < 1 eV, and it only incurs a relative error of up 
to a few percent when 1 < \p\ < 2 eV. 

At this point, it is convenient to define the poten- 
tial Vq = — jti/e, where e > is the proton charge, 
which is associated with the quantum- mechanical effects 
of graphene's band structure [23], and relate it to the in- 
duced charge density per unit area on doped graphene, 
<7 = —en, via the Eq. {T]), 

[dilog (1 + e ~^ v Q) - dilog (1 + e 0eV <*)] ,(2) 



* (kvfPY 

where dilog is the standard dilogarithm function [23l ]. 
One can finally use the definition of differential capac- 
itance per unit area, Cq = da/dVQ, to obtain from Eq. 
([2]) the quantum capacitance of a single layer of graphene 
as [II] 



C Q = — ln[2cosh(/?eVQ/2)] 



(3) 



where we have defined the characteristic length scale for 
graphene, 



(4) 



with the value of X Q w 18 nm at RT. Note from Eq. (JSJ) 
that graphene's quantum capacitance grows practically 
linearly with Vq when this potential exceeds the thermal 
potential, Vth = l/(e/3), having the value of ks 26 mV at 
RT. 

We further assume that an upper surface of graphene 
is exposed to a thick layer of a symmetric z:z electrolyte 
containing the bulk number density per unit volume, N, 
for each kind of dissolved salt ions. Taking advantage 
of planar symmetry, we place an x axis perpendicular to 
graphene and pointing into the electrolyte. The theory 
developed by Borukhov et al. [l9l. |20T| to model finite ion 
size uses the mPB equation for the electrostatic potential 
V(x) in the electrolyte at a distance x from graphene, 
given by 



d 2 V 
dx 2 



Air 



zeN 2 sinh ((3zeV) 
e l + 27sinh 2 (/3zeV/2)' 



(5) 



where z(— 1) is the valency of ions, e is relative dielec- 
tric constant of water (s=s 80, assumed to be constant 
throughout the electrolyte), and 7 = 2a 3 N is the pack- 
ing parameter of the solvated ions, which are assumed 
to have same effective size, equal to a [H, [2(|. We note 
that the standard PB model is recovered from Eq. ([5]) in 
the limit 7^0 fliij ]. By assuming the boundary condi- 
tion V(x) = (and hence dV/dx = 0) at x — > 00, deep 
into the electrolyte bulk, Eq. can be integrated once 
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giving a relation between the electric field and the po- 
tential at a distance x from graphene. Assuming that 
graphene is placed at x = 0, one can use the boundary 
condition at the distance d of closest approach for ions 
in the electrolyte to graphene, 



ax e 



(6) 



to establish a connection between the induced charge 
density on graphene, a, and the potential drop, Vd = 
V(d), across the EDL as 



The total potential, Va, applied between the reference 
electrode in the electrolyte and graphene can be written 
as 



V A = V pzc + V ci + V D + Vq, 



(8) 



where V™ 



(Wg r — VFrcf) /e is the potential of zero 



charge [13l ] that stems from difference between the work 
functions of graphene and the reference electrode, W gl - 
and Wrcf, respectively, and V c \ = Aitha/e' is the poten- 
tial drop across a charge free region between the compact 
layer of the electrolyte ions condensed on the graphene 
surface, having the thickness h on the order of the dis- 
tance of closest approach d [13, [HI, and with e' < e 
taking into account a reduction of the dielectric constant 



of water close to a charged wall [25j. In our calculations, 
we shall neglect these two contributions to the applied 
potential in Eq. ([8]) because V pzc merely shifts the zero 
of that potential, while a proper modeling of V c \ involves 
large uncertainty [20J, 124J . However, usually the effects of 
V c \ can be considered either small [201 ] or incorporated in 
the mPB model via saturation of the ion density at the 
electrolyte-graphene interface for high potential values 
[l9| . Consequently, Vd and Vq represent the two main 
contributions in Eq. (|8|), with the Vd being the surface 
potential of graphene that shifts its Dirac point, and Vq 
being responsible for controlling the doping of graphene 
by changing its chemical potential. Finally, we note that 
all results of our calculations will be symmetrical relative 
to the change in sign of the applied potential because of 
our assumption that the effective sizes of the positive and 
negative ions are equal [l9L |2C| , but this constraint can 
be lifted by a relatively simple amendment to the mPB 
model 

Using the relation Va — Vq + Vd we obtain the total 
differential capacitance of the electrolytically top gated 
graphene as 



cr 



dV A 
da 



dV Q , dV D 



da 



da 



C? 



c 



D 



(9) 



where Cq(Vq) is given in Eq. and Cd(Vd) = da/dVo 
is the differential capacitance per unit area of the EDL, 
which can be obtained from Eq. (0) as 0, [24| , 



Cr 



sinh ((3ze\VD\) 



4ttA 



D 



[1 + 2 7 sinh 2 (/3zeV D /2)] «/| ln[l + 2 7 sinh 2 ((3zeV D /2)] 



(10) 



with Ad = \/e/(8TT0z 2 e 2 N) being the Debye length of 
the EDL [l3[- Note that, in the limit of a very low 
potential Vd, and hence for low density of ions at the 
graphene-EDL interface, one can set 7 — * in Eq. (fT0|) 
to recover an expression for the EDL capacitance arising 
in the standard PB model [HI], 

C D m -r-V- cosh (pzeVo/2) . (11) 

47T AD 

We further note that, while Eq. (fTTj) implies an un- 
bounded growth of the EDL capacitance with Vd in the 
PB model, Eq. (JTUJ) suggests a non-monotonous behavior 
that will eventually give rise to a saturation of the total 
gate capacitance at high applied voltages. 

III. RESULTS 

Given the vast ranges of various parameters in our 
model, it is of interest to define reduced quantities. With 



the thermal potential Vth = 1/ (e/3) , all potentials can 
be written as V = V/V t h- While typical regimes of 
graphene doping require only \Vq\ ^$ 50, we shall extend 
this range in our calculations up to about \Vq\ « 100 
to represent doping levels attained in recent experiments 
in electrolytic environment Q. Next, referring to Eq. 
([2]) , we define the characteristic number density of doped 
charge carriers in graphene by no = (2/ir)/ ((3Tivp) , 
which has the value of no ~ 10 11 cm~ 2 at RT. There- 
fore, defining the reduced density by n = n/no, and 
hence a = a/(eno) = n, we note that |n| may reach 
up to around 10 3 0, [l2T |. It is worthwhile mentioning 
that graphene's characteristic parameters Xq and no are 
related via eA^Agno = 1, where Xb = (3e 2 /e is the Bjer- 
rum length of the aqueous environment, taking the value 
of Xb ~ 0.7 nm at RT [ll|. Furthermore, it follows from 
Eq. that the natural unit of capacitance for this sys- 
tem is Cq = eno/Uth = A^ 1 , taking the value of Co ~ 0.6 
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FIG. 1: Reduced density n of doped charge carriers in 
graphene versus the reduced applied voltage Va for several 
values of the reduced salt concentration N in a. NaF aqueous 
solution. 



FIG. 2: Reduced density log 10 (n) of doped charge carriers in 
graphene versus the reduced salt concentration log 10 (JV) for 
several values of the applied voltage Va (in Volts) for a NaF 
aqueous solution. 



/iF/cm 2 at RT. Turning now to Eq. (0), one can define 
the characteristic number density of ions per unit volume 
in the solution by 



vr 2 2e 2 
N ° 2 U ° Xb „ep(hv F r 



(12) 



which takes the value of N Q w 1.08 x 1CT 6 nm~ 3 w 1.8 
/iM at RT. Defining the reduced concentration of ions in 
the bulk of the electrolyte by N — N/Nq, it would be 
of interest to explore a broad range of its values, e.g., 
1CP 1 < N < 10 5 . Finally, in order to estimate the pack- 
ing parameter, we write 7 = vN and take a = Xb to 
obtain v = 2\%No f» 7 x 10 -7 . With these definitions, 
Eqs. ([2|) and([7]) now read, respectively, 



dilog ( 



1 + e 



-Vr 



dilog ^1 + c 



Vq 



(13) 



2;/ 



In 



+ 2WVsinh 2 (zVd/2) sgn(Vb). (14) 



We now use Eqs. (|13|) and (fl4|) in conjunction with the 
relation Va — Vq + Vd to eliminate the potential compo- 
nents Vq and Vd , and to evaluate the reduced density of 



doped charge carriers in graphene, n = a, as a function 
of the reduced applied voltage Va and the reduced salt 
concentration N. The results are shown in Figs. 1 and 
2, covering the following ranges: |n| < 10 3 (correspond- 
ing to |n| < 10 14 cm" 2 ), \Va\ < 60 (corresponding to 
\V A \ ^ 1.6 V), and 10" 1 < N < 10 5 (corresponding to 
0.18 fjM < N < 0.18 M). In Fig. 1 one notices a strong 
dependence of n on the applied potential for Va greater 
than about 30, which gives approximately equal rates of 
change for each salt concentration at the highest values 
of the applied potential. On the other hand, at the lower 
applied potential values, there exists a much stronger de- 
pendence on the salt concentration, which is revealed in 
Fig. 2, showing log 10 n versus log 10 TV for several applied 
voltages. Indeed, one notices a very strong sensitivity 
of the doped charge carrier density in graphene to the 
salt concentration for applied voltages Va ^ 0.4 V in the 
range of salt concentrations N <■ 1 mM. Even though this 
sensitivity seems to be the strongest at the lowest applied 
voltages, one should bear in mind that the electrical con- 
ductivity in graphene becomes rather uncertain around 
its minimum value, which extends up to doping densi- 
ties about n w 10 11 cm -2 [2a, [23. Therefore, it seems 
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that 0.1 ^ ~ 0.4 V would be an optimal range of 
applied voltages for sensor applications of the electrolyt- 
ically top gated graphene in probing salt concentrations 
in the sub-microMole range. 

Next, moving to the capacitance of electrolytically top 
gated graphene, C, we note that the reduced capaci- 



tances, Cq = Cq/Cq and Cd — Cd/Cq, are obtained 
from Eqs. © and (TO]) as 



Co = 2 In 



2 cosh 



(Vq/2) 



(15) 



N sinh 



1 + 2^iVsinli (zV D 



>/2) 



(z\Vd\) 



In 



1 + 2vN sinh" 



(zVd/2) 



(16) 



showing that Cq and Cd are comparable in magnitude 
for vanishing potentials when the salt concentration is 
N ~ 1. Moreover, referring to Eq. ([9]) as an electrical 
model where graphene and the EDL act as a series con- 
nection of capacitors, it follows that graphene's quantum 
capacitance Cq will be promoted as the dominant con- 
tribution to the total gate capacitance as the salt con- 
centration increases. 

We now use the equality of the right-hand-sides in Eqs. 
([T3|> and QUI along with the relation V A = Vq + V D to 
eliminate Vq and Vd, and to evaluate the reduced quan- 
tum capacitance of graphene from Eq. (|15| . as well as 
the reduced capacitance of the EDL from Eq. (fill)) as 
functions of the reduced applied voltage Va- Results are 
shown in Fig. 3 along with the total reduced capacitance 
of the system based on Eq. for two reduced salt con- 
centrations, N = 1 and 10 5 (corresponding to N ~ 1.8 
fjM and 0.18 M, respectively). We show our results for 
the reduced applied voltages up to Va ~ 120 in order to 
elucidate the effect of saturation in the total capacitance 
that occurs at Va « 85 V (corresponding to Va ~ 2.21 
V) for N = 1 and at Va ~ 75 (corresponding to Va ~ 
1.95 V) for N = 10 5 . As can be seen from dotted curves 
in Fig. 3, showing a non-monotonous dependence of the 
EDL capacitance on the applied voltage, the saturation 
effect in the total capacitance of the electrolytically top 
gated graphene is a consequence of the steric effect of the 
electrolyte ions that are crowded at the graphene surface 
at high applied voltages [20]. Even though the voltages 
where the saturation occurs are relatively high, they may 
still be accessible in experiments on graphene. Further- 
more, we see in Fig. 3 that, at intermediate applied volt- 
ages, the rate of change of the total capacitance follows 
closely that of the quantum capacitance, with the value 
~ 23 /jF/ (Vcm 2 ) that is commensurate with recent mea- 
surement [ll|. At the lowest applied voltages, one notices 
in Fig. 3 a "rounding" of the total capacitance as a func- 
tion of voltage for low salt concentrations, which comes 
from the EDL capacitance. Such rounding is observed in 
the recent experiment [Til ]. 

As a consequence of the vast differences between the 




FIG. 3: The dependence on the reduced applied volt- 
age Va is shown for: the total reduced capacitance C = 
CqCd/ [Cq + CdJ (solid black lines), graphene's reduced 

quantum capacitance Cq (dashed red lines), and the reduced 
capacitance of the electric double layer Cd (dotted blue lines), 
in the NaF aqueous solutions with the reduced salt concen- 
trations of N = 1 (thick lines) and 10 5 (thin lines). 



capacitances shown in Fig. 3, one expects that there ex- 
ists a broad variation in the way how is the total applied 
voltage Va split between the potential drop Vd across 
the EDL and the voltage Vq pertaining to the change in 
graphene's chemical potential. We therefore display in 
Fig. 4 the variation of the fraction r = Vd /Va as a func- 
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FIG. 4: The dependence of the ratio r = Vd/Va on the re- 
duced salt concentration log 10 (iV) for several values of the 
applied voltage Va (in Volts) for a NaF aqueous solution. 



tion of the reduced salt concentration in the electrolyte 
N for several values of the applied voltage Va- One can 
see that, at low salt concentrations, the potentials Vq 
and Vd are roughly comparable in magnitude, although 
the ratio r increases in favor of the potential drop across 
the electrolyte as the applied voltage increases. How- 
ever, this trend is reversed at high salt concentrations 
and, more importantly, Fig. 4 shows that most of the 
applied voltage is used to increase graphene's chemical 
potential for a full range of applied voltages when salt 
concentration N exceeds, say, mM. Besides its impor- 
tance for applications, the fact that the potential drop 
across the electrolyte remains very small at high applied 
voltages also alleviates concern that a high electric filed in 
the electrolyte may cause the onset of voltage-dependent 
electrochemical reactions on graphene. 



IV. CONCLUDING REMARKS 

We have analyzed the doping of single-layer graphene 
due to application of the gate potential through an aque- 
ous solution of salt using a modified Poisson-Boltzman 
model for electrolyte, and found great sensitivity of the 
induced charge density in graphene to the broad ranges 
of both salt concentration and applied voltage. We have 
further analyzed differential capacitance of the electrolyt- 



ically top gated graphene, and found that its quantum ca- 
pacitance is promoted as the dominant component owing 
to a reduction in the Debye length of the electric double 
layer when the salt concentration increases. In this case, 
very little potential drop appears across the electrolyte, 
and graphene takes most of the voltage drop to shift its 
chemical potential. These findings have several impor- 
tant consequences. 

First, since graphene's conductivity is dependent upon 
its chemical potential via doping density, Eq. |T]), the 
sensitivity of the latter to the salt concentration im- 
plies good prospects for applications in biochemical sen- 
sors, especially for in- vivo electrochemical measurements 
in biological systems owing to graphene's natural bio- 
compatibility. Next, since most of the applied voltage can 
be used to increase the chemical potential of graphene, 
as opposed to a potential drop across the electrolyte, one 
can envision ways to use a very thin top gate (in the form 
of a liquid or solid electrolyte) that requires relatively low 
gate voltage to change the chemical potential (and hence 
conductivity) of graphene in future small scale field ef- 
fect devices with tunable conductivity. Among other as- 
pects of the increased role of graphene's quantum capac- 
itance is reduction of the electrical field across the elec- 
trolyte. This can help reduce the rates of voltage- induced 
electrochemical reactions on graphene's surface, as well 
as improve the mobility of charge carriers in graphene 
by reducing their scattering rates on various impurities. 
Moreover, since quantum capacitance is basically the ca- 
pacitance associated with change in carrier densities in 
graphene, it can be seen as analogous to the junction 
capacitance, and the smaller quantum capacitance could 
in turn lead to faster switching time for graphene based 
devices. 

Many of these advantages of top gating through an 
electrolyte are related to a high bulk dielectric constant 
of the electrolyte, especially in aqueous solutions. So, 
even though the oxide thickness can be reduced down to 
around 2nm in the present generation conventional MOS 
structures, the much higher dielectric constant of water 
in comparison to SiC>2 should provide for a higher gate 
capacitance, translating into much better field effect per- 
formance, as discussed above. However, dielectric con- 
stant of an electrolyte can be significantly reduced close 
to a charged surface 0, [25j], and this issue has yet to 
be discussed in the context of electrolytic top gating of 
carbon nano-structures. 
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